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Abstract 

The low-energy modes of a superfluid atomic Fermi gas at zero temperature are investigated. 
The Bose-Einstein-condensate (BEC) side of the superfluid phase is studied in detail. The atoms 
are assumed to be in only one internal state, so that for a sufficiently diluted gas the pairing 
of fermions can be considered effective in the I = 1 channel only. In agreement with previous 
works on p-wave superfluidity in Eermi systems, it is found that the px + ipy phase represents the 
lowest energy state in both the Bardeen-Cooper-Schrieffer (BCS) and BEC sides. Our calculations 
show that at low energy three branches of collective modes can emerge, with different species of 
dispersion relations; a phonon-like mode, a single-particle-like mode and a gapped mode. A 
comparison with the Bogoliubov excitations of the corresponding spinor Bose condensate is made. 
They reproduce the dispersion relations of the excitation modes of the p-wave superfluid Fermi 
gas to a high accuracy. 
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I. INTRODUCTION 


Over the last several years a sustained interest has been devoted to paired fermion phases 
with unconventional pairing symmetry. Among these the p-wave spin triplet condensate has 
attracted particular attention. The occurrence of p-wave superfluidity in Fermi systems was 
already studied in the sixties by Anderson and Morel [1] in relation to the low temperature 
phase of liquid ^He. In the last decade there has been a renewed interest after the obser¬ 
vation of p-wave Feshbach resonances in ^Li and atoms M, which had raised the 
prospect of realizing p-weive superfluidity in cold Fermi gases. In particular, the possibility 
of controlling the strength of the interatomic interactions via Feshbach resonances has stimu¬ 
lated theory works on the evolution of the superfluidity from the Bardeen-Cooper-Schrieffer 
(BCS) regime to the Bose-Einstein condensation (BEC) of composite bosons [UffTTj. We re¬ 
mark that a p-wave superfluid phase may occur also in the Fermi component of a gaseous 
mixture of ultracold bosons and one-component (spin-polarized) fermions. In that case an 
attractive interaction between fermions can be induced by the exchange of virtual phonons 

mm- 

In this work we explore properties of a three-dimensional spin-polarized Fermi gas in the 
p-wave superfluid phase using a fermion-only model. Since for the p-wave superfluid phase 
the order parameter is given by a three-dimensional object two distinct phases are available 
to the superfluid. A p^. - 1 - ipy phase, which is characterized by the 2 ;-component of the 
pair angular momentum m equal to ±1 and a Px phase with m = 0. At low temperatures 
the system can undergo a phase transition from a Pa, to a p^, -|- ipy superfluid state [7]. 
Moreover, detuning the Fesh- bach resonance the system can exhibit a transition from a 
gapless quasiparticle spectrum for positive values of the chemical potential p (BCS-side) 
to a gapped spectrum for p < 0 (BEC-side) [71 HI]. We will address in particular the 
study of low-energy collective modes in the limit of vanishing temperature. The threshold 
for pair-breaking vanishes in the BCS-side. So the collective modes in general should be 
strongly damped in the BCS-side. For this reason, we will consider collective excitations 
of the Fermi superfluid in the BEC phase only. Collective modes in a p-wave superfluid 
Fermi gas have been already studied in Ref. [12]. However, in that paper calculations 
were performed only for the Px superfluid phase. Moreover, the authors reported dispersion 
relations only for the Nambu-Goldstone mode related to phase-oscillations of the pairing 
held. We aim to do a step further. We consider the system in the px + ipy phase, which 
represents the true ground state, and we extend calculations to collective modes given by 
amplitude huctuations of the pairing held. Furthermore, we derive from the properties of the 
superhuid Fermi gas in the BEC side the structure of the corresponding Bose counterpart 
and the interaction between the composite bosons. Since the order parameter of the fermion 
superhuid phase is represented by a complex vector the boson gas can be considered as a 
Bose spinor condensate. We compare the dispersion relations for the collective modes of the 
Fermi system with the corresponding ones of the Bose condensate. From an experimental 
point of view the observation of the peculiarities of the excitation spectrum may help to 
assess the occurrence of a BEC regime in a p-wave superhuid Fermi gas 1201 . 


II. FORMALISM 

In the imaginary-time functional formalism an ehective action for the pairing held can 
be obtained by using a Stratonovich-Hubbard transformation for the couples of fermion 
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fields (see, e.g., Ref. [ 21 ]). After integration of the fermionic part the effective action in 
momentum representation is given by 


S.„(A,A*)= -t ^ /'’dr[A-(P,k,T)l/-‘(k,k')A(P,k',T) 

-5i:^-^Tr(e„Ar, 


( 1 ) 


where is a 2 x 2 diagonal matrix, whose elements are the imaginary-time propagators for 
independent fermions with chemical potential /i, and the matrix A 


A(P,k,r) 


0 A(P,k,r)\ 

AAP,-k,r) 0 ) 


represents a complex bosonic held, which is periodic in the imaginary time interval (0, (3 = 
1/T) (units such that h = ks = I are used). 

In Eq.(|^ the matrix l/“^(k, k') is the inverse of the fermion-fermion interaction 

< P/2 - k, P/2 + k\V\P/2 + k', P/2 - k') >= ^Vi\k - k'l) (2) 

where k = (ki — k 2)/2 and k' = (ks — k 4)/2 are the relative momenta of a pair of fermions, 
and P = ki+k 2 = k 3 +k 4 is its center-of-mass momentum. In the above equation l/(|k—k'|) 
represents the space-Fourier transform of the interaction and V is the normalization volume. 

The pairing held at equilibrium, Ao(P, k, r), is evaluated within the saddle-point approx¬ 
imation to the ehective action S'e//(A, A*), while for the huctuations of A(P, k, r) about its 
equilibrium values we adopt a gaussian approximation. The equation for Ao(P,k, r) reads 

Ao(P, k, r) = ^ R(|k - k'|)^(i^)(k' + P/ 2 , k' - P/ 2 , r - r' = 0), (3) 

^ k' 


where + P/2,k' — P/2,r — r' = 0) is the equal-time anomalous propagator for 

fermions interacting with the pairing held A |22j . 

The solutions of Eq. (|^ with P 7 ^ 0, correspond to a breaking of the space-translational 
symmetry, i.e. the LOFF phase [23]. Here, only solutions with vanishing center of mass 
momentum are considered. In this case the gap equation becomes 


Ao(k) 


1 

V 


Erdk-k'D/i: 

k' P n 


Ao(kO 

+ ik' + ‘Ao(k')Ao(k') 


( 4 ) 


Now, we expand the ehective action Seff{A, A*) up to the second order in the huctuations 
of the pairing held 

cr(P, k, r) = A(P, k, r) - Ao(k) 
and for the ehective action of the huctuations we get 


6Seffia,a*) =-j Y. [ dr[a^(P,k,r)R ^(k,k')a(P,k',' 

^ P,k,k' 

“tH / drdr'[a^(P,k,r)A(r-r',P,k)a(P,k,r') 

4p,k^o ^ 


( 5 ) 
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with the vector a(P, k, r) and the matrix V ^(k, k') given by 


d(P,k,r) 


a*(P,k,r) ) 


and 


t>-'(k,k') 


P-Hk,k') 0 \ 

0 P-Hk,k') ) • 


In Eq.([^ we have introdnced the 2x2 matrix 
A{t — t', P, k) whose elements are 


- r', P, k) = - t\ P/2 + - r, P/2 - k) , 

dl 2 , 2 (r - r', P, k) = - r', P/2 + \i)g^'^^\T' - r, P/2 - k), 

2 li, 2 (r - r', P/, k) = g^^^\T - r', P/2 + \i)g^^^\T - r', P/2 - k), 

A 2 ,i{t - t', P, k) = - r', P/2 + - r', P/2 - k), (6) 

where the Green fnnctions in the r.h.s. are calcnlated for the gas at eqnilibrinm. Then, 
they are diagonal in the momenta. From general properties of the Green fnnctions [22] the 
following symmetry relations ensue 


A 2 , 2 {r - r, P, k) = Ai_i(r' - r, P, k) 


and 


212 ,i(r - r', P, k) = AlAr - r, -P, -k). 

Equation ([^ shows that the propagator for the bosonic held [cr*(P, k, r), (j(P, k, r)]l 
obeys the integral equation 


Dir - r', P, k, k') = -t/(|k - k'|)5(r - r') 


/p dri[P(|k-k”|)A(r-ri,P,k”)T)(ri-r',P,k”,k')] . (7) 


In order to calculate the dispersion relations of collective modes the above equations 
should analytically be continued to real times. This can simply be accomplished by substi¬ 
tuting in the expressions for the matrix elements Aij the imaginary-time Green functions 
with their real-time counterparts . Then, we recast Eq.([^ in the real-frequency represen¬ 
tation 

B(w. P. k, k') = I{^(k, k') - / i: '^(k. k” P, k" )B(w. P, k”, k'). ( 8 ) 

y V 

Here, we consider a gas in the limit of vanishing temperature T = 0. In this case the 
matrix elements P, k) are explicitly given by 


and 


Hi,i(a;,P,k) 


n^(k)M^(k) n^(k)ni(k) 

u -|- E+(k) -|- E_(k) — ie u — E+(k) — E_(k) -|- ie 


A+(k)A_(k) E+(k) + E_(k) 

2 E+(k)E_(k) a ;2 - (E+(k) + E_(k) + ief ’ 


(9) 

( 10 ) 
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while the remaining matrix elements are determined by the symmetry relations 


A2,2{uj, P, k) = P, k), 

A2,i{^, P, k) = -P, -k). 

In the above equations we have used the notations 

u±(k), u±(k), A±(k), E±(k) = u(P/2 ± k), u(P/2 ± k), Ao(P/2 ± k), E{P/2 ± k), 


where i?(P/2 ± k) = /2 E'k) + |Ao(P/2 ± k)P is the quasiparticle energy, 


and 


u^{P/2±\i) 


u2(p/2±k) 


1. e(P/2±k) s 

2V ^E(P/2±k)^ 

1 . i{P/2 ± k) ^ 

2V E(P/2±k)^ 


with ^(P/2 ± k) = (P/2 ± k)V2m - /i. 

The dispersion relations of elementary excitations are given by the real poles of the prop¬ 
agator, Eq.(j^, for given values of P. In the vicinity of the poles the finite term 2I/(|k, k'|)/V 
in Eq. ([^ can be neglected and the poles are determined by equating to zero the determinant 
associated to the homogeneous set of equations 


D{u, P,k, k') 


1 

V 


X:P(k,k”)A(u;,P,k”):D(u;,P,k”,k'). 


The above equations admit solutions of the form 

D{u, P, k, k') = P, k) X Xb{(^, P, k') 


with the two-dimensional vector x(a;, P, k) = Xa{^^, P, k) obeying the equations 

P, k) = -/ E P, k')x(wP, k'). (11) 

In the present work we are mainly interested in the global features of the collective modes. 
Then, a simplified interaction between fermions may be satisfactory for our purpose. Firstly 
we assume that in the expansion of the interaction in Legendre polynomials 

k(|k - k'l) = E 

the term with I = 1 gives the dominant contribution to the pairing of fermions. We observe 
that, because of the Pauli principle, only the components with odd values of the relative 
angular momentum I can be effective. In addition, we adopt for Vi{k,k') the schematic 
form Vi{k,k') = —Xkk', with A > 0, for k,k' < k^., and Vi{k,k') = 0 otherwise. It is 
convenient to express the interaction strength A and the momentum cut-off k^ in terms of 
the scattering volume Vi and the second coefficient, ri, in the effective-range expansion of 
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the scattering amplitude [2U |25] . By exploiting the low energy expansion of the scattering 
matrix [9l [15], [26] one can obtain the relations 



^ _ _47r^^ _ _4 


( 12 ) 


where = k‘^/2m. The diluteness condition nrf << 1 implies [kc/kp) » 1, where n is 
the density of atoms and kp \s the Fermi momentum. 

The equation (|^ in the limit T = 0 becomes 


Ao(k) 





dk! 

(27r)3 


k'Y^ii!) 


Ao(kO 

2E(k') ’ 


(13) 


and, in order to determine explicity the held Ao(k) together with the chemical potential /i 
the equation hxing the fermion density has to be added 


Up 


^k 1 / m \ 

{2nf2\ E(k); ’ 


(14) 


The energy per fermion in the superhuid phase is given by the expression 


Ep = — 


dk 


npJ {2Tif2 




2E(k) 


+ h) 


(16) 


which coincides with the usual expression of the BCS theory [22], apart from a factor 1/2 
due to the absence of degeneracy for the Fermi gas in the present case. 


Equation (13) suggests solutions of the form Ao(k) = Aq ■ k 6 '(/cc — k), where the vector 
Ao can be of two different species a real vector (pa,-phase) or a complex vector 

{px+ipy-ph.a.se). In the latter case Aq can be decomposed according to Aq = Ai+iA 2 where 
Ai and A 2 are two real vectors of equal magnitude and each other orthogonal. Moreover, 
vectors obtained by global gauge transformations and/or by rotations, in the orbital space, 
of a particular choice of Aq represent equivalent solutions for the equilibrium pairing held. 
Explicit calculations show that the Px + ipy phase is the lowest-energy state [U [71I2] ■ 

With the adopted interaction between fermions, Eq. becomes 


x(L^,p,k) = \^l-'£Yr(kk 


dk' 


i'r”(k')-4(w,P,k')x(w,P,k') 


The amplitudes x(a;, P, k) can be conveniently expanded as 

x(a;,P,k) = ^x(-)(a;,P)A;Fr(k), 

m 

and we obtain for the m-components the set of coupled equations 


x^’"Hu,,P) = '£ 


Ai/y I -^^k’X’"(k)A(c,,p,k)Yr'’{k) 


X 


(m') 


U,,P). 


(16) 
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III. RESULTS 


In this paper we are concerned with the fluctuations of the pairing held in the stable 
phase Px + ipy The ground-state energy and the chemical potential are determined by the 
magnitude Aq = -^/Aq ■ Aq alone. In order to calculate Aq we choose a particular direction 
for the vectors Ai and A 2 , say Ai = (Ao/\/2, 0,0) and A 2 = (0, Ao/\/2, 0)_. Moreover, we 
observe that once the value of the cut-off kc is fixed, the scaled quantity Aq = Aokp/ep, 
where ep = k"^/2m is the Fermi energy, depends only on the dimensionless effective coupling 
constant 

A = 2mXkp , 

which is a function of the product vikp and of the scaled cut-off kc/kf. In terms of the 
parameter A the equation for Aq reads 


1 



dk 1~2 . 2rn \ — k) 

- — --k^sin^ieA X — , , ^ , J , 

(2^)" 2 2^i\k) + lAfk^stn^e^) 


(17) 


with ^{k) = ^{k)/ep and momenta expressed in units of the Fermi momentum: k,kc = 
k/kp, kc/kp. 

The dispersion relations for collective modes, n;(P), are given by the frequencies, for 
which the determinant of the matrix 


1 - At{u:, P) = P) 


with 


A. 


T(iJ) 



k%^(k)Aju,p,k)Yr'*(k) 


(18) 


vanishes for given values of the momentum P. Similarly to the_ energy gap the scaled 
frequency uj(P)/€p depends only on the dimensionless parameter A and the scaled cut-off 
kc 

With the particular choice for the plane where the vector Aq lies and for the phase of 
its components, the spatial 50(3) symmetry and the gauge U{1) symmetry are broken. We 
notice that in the present case the spins of the fermions do not play any role since they are 
frozen along a fixed direction and there is no coupling between spin and orbital degrees of 
freedom. As a consequence, we expect the occurrence of Nambu-Goldstone modes related 
both to the phase fluctuations of the m = 1 component of the pairing held and to the 
fluctuations of the amplitudes of the pairing held. Actually, by exploiting the gap equation 


one can see that Eqs. (16) show solutions with a; = 0 for vanishing momentum P. However, 


this does not give rise in general to undamped phonon-like excitations. This occurs because 
for particular directions of the momentum the energy gap vanishes, so that in the BCS 
side the phonon-like modes can merge in the continuum of ungapped two-quasiparticle 
excitations. 

The threshold for the quasiparticle continuum is given by the minimum value of the 
sum E(P/2 + k)) + E(P/2 — k). The most favorable case is when P is ’’parallel” to the 
complex vector Aq. In this case the threshold is given by AoP/a/ 2, and only phonon-like 
modes with a phase velocity such that vs/vp < Ao/( 2 )T 2 ^ where vp is the Fermi velocity, 
can propagate without damping. Explicit calculations show that this requirement can be 
satisfied only for large energy gaps. The corresponding values of the interaction strength are 
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l/(k’pV,) 


FIG. 1: The scaled value of the p-wave energy gap as a function of l/{vikp) for two values of the 
momentum cut-off: kc/kp = 10 (red circles) and kc/kp = 20 (blue circles). 


in the region close to the BCS-BEC transition. We should point out that our calculations 
are essentially performed within the framework of a mean-held theory. The extension of 
the present approach to a critical region, where huctuations of the pairing held could be 
important, can be rather questionable. Therefore, we do not discuss collective modes in the 
BCS side and focus our attention on the collective excitations of the Fermi gas inside the 
BEG phase. 

The scaled energy gap Aq and the chemical potential together with the energy per particle 
are reported in Figs. (1) and (2) respectively as functions of the parameter l/{vikp) for two 
diherent values of the cut-oh, kc = 10, 20. The curves for kc = 10 show a steeper behavior. 
This is simply due to the fact that the ehective coupling constant A is larger for smaller kc 
at a given value of l/{vikp). Figure (2) shows a remarkable closeness between the chemical 
potential and the energy per particle, the two sets of points are practically indistinguishable. 
This indicates that the pairs of fermions behave as a weakly interacting Bose condensate 

m- 

Now we turn to the dispersion relations of collective modes. In order to bring out any 
anisotropy in the propagation of the pairing-held huctuations we consider two directions of 
the wave vector, orthogonal each other: P directed along the z axis, i.e. perpendicular to 
the complex vector Aq, and P lying in the (x, y) plane, parallel to the x axis for dehniteness. 

In the hrst case, the quasiparticle energies E(P/2 ± k) do not depend on the azimuthal 
angle, 0, of k. Then, the diagonal matrix elements Aj j(a;,P,k) are independent of 0 as 
well. Whereas the oh-diagonal elements P, k) and ^2,1(0;, P, k) are multiplied by the 

factors and respectively. As a consequence, the amplitude with m = 0 is decoupled 
from the remaining amplitudes. In addition, x^^^(a;,P) is coupled only to xi 
the amplitudes Xi and X 2 ^^(a;,P) separately obey single equations. In the other 

case, where P is parallel to the x axis, the situation is quite different. The amplitude with 
m = 0 is still independent of the remaining ones. Whereas, for the amplitudes with m = ±1 
we have to solve four coupled equations. However, explicit calculations show that the matrix 
elements, which couple x-”^^(ci;,P) to X'~™'^(a;,P) and Xi”*^(a;,P) to X2”^^(a;, P), are about 





FIG. 2: Energy per particle (circles) and chemical potential (dashed lines), in units of the Fermi 
energy, as a function of l/{vik'p) for two values of the momentum cut-off: kc/kp = 10 (red color) 
and kc/kp = 20 (blue color). 


three orders of magnitude smaller than the other ones. Practically, the structure of the set 
of equations is the same in both the cases, but the values of the coefficients are different in 
general. 

For the effective coupling constant we £x the value corresponding to vikp = 0.1 and 
consider two different values of the scaled cut-off kc = 20 and kc = 10. The values, in units 
of ep, of the energy gap, chemical potential and energy per particle are 

Ao = 0.484, /i = -1.70, 

for kc = 20 and 

Ao = 0.821, /i = -4.28, 

for kc = 10. 

In Figs. (3) and (4) the dispersion relations a;(P) are reported for the two cases considered 
and for modes propagating along a direction perpendicular to the equilibrium pairing held 
(z axis). For modes propagating along the x axis the dispersion relations are very similar, 
the sets of points are indistinguishable from those of Figs. (3) and (4). Practically the 
propagation of the collective modes can be considered isotropic, the difference between the 
phase velocities for two orthogonal directions amounts to less than 1% on average. 

Figs. (3) and (4) show some remarkable features of the spectra of collective excitations. 
In addition to the usual ungapped phonon-like branch related to phase huctuations of the 
equilibrium pairing held (m = 1), there are a gapped spectrum for the m = — 1 amplitude 
and a parabolic dispersion relation for the m = 0 amplitude. From these properties of 
the dispersion relations we can gain some useful insight on the structure of the BEG and 
the interaction between composite bosons. The wave-function of the relative motion of the 
couple of fermions in the BEG side is an eigen-function of the relative angular momentum 
with I = 1 and m = 0,±1 [19]. If the details of the internal structure are neglected 


E 

N 


— = -1.72 


E 

N 


— = -4.34 
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FIG. 3: Dispersion relations of collective modes of the superfluid Fermi gas (circles) and excitations 
energies of the corresponding Bose condensate (dashed lines) as a function of the scaled wave vector 
for l/{k\vi) = 10. The scaled momentum cut-off is kc/kp = 20. Red, blue and green colors 
correspond to m = 1, —1,0 oscillations respectively (see text). 



FIG. 4: Same as Fig. 3, but for kc/kp = 10. 


the composite bosons can be considered point-like particles with spin S = 1. We have 
studied their equilibrium state and their dynamics by means of the time-independent and 
time-dependent Gross-Pitaevskii equations (see, e.g., Ref. 1251 ). Since the spectrum of 
elementary excitations is isotropic in orbital space, it is not necessary to introduce any 
coupling between spin and orbital degrees of freedom. The simple interaction 

< qi,mi;q2,m2|R|q3,m3; q4,m4 >= 5qi+q2,q3+q4(7o+ 7i < mi,m2|Si ■ 821 ^ 3,7714 >) (19) 

has been used. For an extensive review of spinor Bose-condensates see, e.g., Ref. [28] . 

Since the Fermi system under consideration is in a + ipy superfluid phase, we should 
expect that the corresponding Bose held at equilibrium, contains only the 777 = 1 

component ( ferromagnetic phase ), which can be assumed real and positive. In this phase, 
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from the time-independent Gross-Pitaevskii equation one obtains for the difference between 
the chemical potential and energy per particle of the composite bosons the relation 

= + ( 20 ) 

with hb = 2 p, E/Nb = 2E/N and = hb = n/2. The requirement that represents 
a state of stable equilibrium, implies that the inequalities 71 < 0 and 70 + 7 i > 0 should be 
fulhlled [ 2 H]- 

As we expect a weak interaction between composite bosons, the Bogoliubov approxima¬ 
tion to the solutions of the time-dependent Gross-Pitaevskij may be satisfactory. Within 
this approximation one obtains the following dispersion relations for the eigen-modes of the 
spinor Bose condensate [28]: a Nambu-Goldstone mode 


4'A) = AhAMmhT+TM, 


a single-particle-like mode 


Ef(P) = iB{P) 


and a gapped mode 


where eB{P) = P‘^/{2M, 


B, 


E^-^\P) = eB{P)-2^,nB, 
with Mb = 2m is the boson kinetic energy. 


The scaled energies E^^\P)/ep depend on the coupling constants only through the di¬ 
mensionless quantities P* = 'yiUB/^F, with z = 0,1. The values of P^ can be determined from 


Eq. (20) and from the gap observed in the dispersion relation for the m = —1 fluctuations. 
Figs. (3) and (4). They are 


To = 0.135, Pi = -0.045 

for kc = 20 and 

Po = 0.349, Pi = -0.116 

for kc = 10 . 

In Figs. (3) and (4) the dispersion relations for the eigen-modes of the Bose condensate 
have been added. We observe that they reproduce the corresponding excitation spectra of 
the superfluid Fermi gas with noticeable accuracy. The phase-velocity_ of the phonon-like 
mode, in units of the Fermi velocity, is 0.150 for k^. = 20 and 0.241 for k^ = 10. 

A remark is required about the symmetry properties of both the Fermi superfluid and the 
Bose spinor-condensate counterpart. The superfluid px + ipy phase and the corresponding 
ferromagnetic phase for the BEG break both the rotational symmetry and the U{1) sym¬ 
metry under global gauge transformations. However, U{1) transformations can make up for 
spatial rotations about the original z axis, i.e. the ground state displays a ’’spin-gauge” 
symmetry [291 - I3T] . So that a U{1) symmetry is still maintained. This reduces to 3 the 
number of spontaneously broken group-generators. As a consequence 3 Nambu-Goldstone 
modes should occur. Our calculations show a phonon-like mode and a single-particle-like 
mode, in addition to a gapped mode. According to Ref. [32] the hrst mode should be 
counted once, whereas the second should be counted twice. Then, the Goldstone theorem is 
still fulhlled. 
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IV. SUMMARY AND CONCLUSIONS 


For a sufficiently diluted atomic Fermi gas with only one component, the interaction 
between atoms can be considered effective only in the p-channel. The transition from the 
BCS to the BEC regime can be obtained by tuning the strength of the interatomic interaction 
via a Feshbach resonance. In our approach to assess the occurrence of a superffuid phase in 
the Fermi system and to determine its properties, only fermionic degrees of freedom come 
into play. The lowest energy state is a + ipy superffuid phase in both the BCS and BEC 
sides. This is in agreement with ref. [1] for the BCS regime and with ref. [7] for the BEC 
regime. In the latter paper calculations were performed within a coupled fermion-boson 
model for a p-wave Feshbach resonance. 

We expect that the collective modes are strongly damped in the BCS side. Thus, we 
have studied the low-energy excitation spectrum only in the BEC phase. In this phase 
the dispersion relations of the collective modes show peculiar features, which might also be 
relevant for the transport properties of the system. Besides an excitation mode with the usual 
phonon-like dispersion relation, a gapped mode and a single-particle-like mode can occur. 
However, these features are consistent with the Goldstone theorem. Moreover, when the 
superffuid Fermi gas is well inside the BEC phase the chemical potential and the energy per 
particle approach each-other. This suggests that the pairs of fermions behave as a Bose- 
Einstein condensate with a weak repulsive interaction acting between composite bosons. 
Finally, a comparison has been performed between the dispersion relations of the Fermi- 
gas collective modes and the corresponding ones of a Bose spinor condensate calculated by 
means of the time-dependent Gross-Pitaevskii equations. They coincide to a high accuracy. 
Then, there is a tight similarity between a p-wave superffuid Fermi gas in the BEC side and 
a Bose spinor condensate for both the equilibrium properties and the dynamical behavior 
at low energy. 

A final comment is in order. In a mean-field approximation retardation effects, arising 
from the renormalization of the effective interaction between fermions in a medium, are 
neglected. In Ref. [22] it has been shown that the energy and momentum dependence of the 
interaction plays a significant role in determining the p-wave pairing gaps in a superffuid 
Fermi gas. Then, we should expect that, from a quantitative point of view, the inclusion of 
such effects might affect the results, presented in the present paper, appreciably. 
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